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MULTIPLE EXOSPHERIC TEMPERATURES FOR VENUS AND MARS 


ABSTRACT 

Solutions of the heat transport equations for each of the important 
neutral species in the upper atmosphere of Venus and Mars have shown that 
multiple neutral temperatures must be considered to properly understand 
tiie distribution of atomic hydrogen. Temperature differences between CO^ 
and the light atomic species, H, D, and He in excess of 100 °K for Venus 
and 50 °K for Mars are obtained above the exobase. The magnitude and 
separation of these temperatures depend strongly on the relative proportion 
of light neutrals to carbon dioxide in the lower atmosphere and on the degree 
of turbulent mixing in the lower thermosphere (above 100 km) . Since the 
atomic hydrogen effusion flux is strongly dependent on the exobase tempera- 
ture, the exospheric hydrogen density is found to be very different from that 
based on the C0 2 temperature. The exospheric atomic hydrogen density is 
found to have a power law dependence on the proportion of neutral helium 
present in the atmosphere. 


MULTIPLE EXOSPHERIC TEMPERATURES FOR VENUS AND MARS 


The detailed knowledge we have of the neutral gas temperature above 
100 km lias been obtained from the Mariner and Venera series of flights to 
Venus and Mars. The Lyman a measurements made by Barth et. al (1968 , 1971) 
and Kurt, et. al (1968) when coupled with an appropriate model of an exosphere, 
Wallace (1969) can be made to yield hydrogen temperatures and densities at 
high altitudes. Additional data can be obtained from electron density profiles 
such as the familiar data obtained by Mariner V, Fjeldbo and Eshleman (1969). 

The neutral gas temperature can be inferred from the slope of the density pro- 
file above the ionization peak, Figure 1, if the atmospheric composition or 
mean mass is known. The procedure is valid because below about 180 km altitude 
the electron, ion, and neutral gases exhibit strong collisional coupling. This 
coupling both forces a common temperature T n - T^ and minimizes the effects of 
ion-electron diffusion on the electron density scale height. Thus, the small 
scale height above the ionization peak has been used to infer a mainly C0 9 + 
lower ionosphere in a CC^ atmosphere with an exospheric temperature of 650 °K 
for Venus and 400 °K for Mars. Since the rapid change in scale height above 
200 km on Venus clearly indicates a change in composition from C0 ? to light 
neutral gases, do the exospheric temperatures obtained from the electron density 
data also apply to the light atomic species He, H, and D at high altitudes? 

It will be shown that the general answer is no; the temperatures for CO., and 
the lignt gases are different except for certain special cases of atmospheric 
composition. 

The computation of the temperatures was performed by the methods given 
in a recent paper by et. al (197V with the addition of an independent 
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heat flow equation for the light neutral gases, H, D, and He. Since the 
thermal behavior of these gases is quite similar it was decided to assign 
a common temperature to H, D, and lie, but not to the heavier molecular gases 
such as CO^. For example the thermal conductivities of H, D, and He arc 

0.73 0.73 0.75 _i _ 1 .1 

16. 4T , 11.7T , and 20. 9T erg cm sec °K respectively, and 

1-23 

for C0-, it is 1.50T , Tsederberg (1965). The heat input and loss mechan- 

isms for H, D, and He are quite similar since they are dependent on electron- 
neutral, ion-neutral and neutral -neutral collisions rather than radiative 
processes as with CO^. Thus, even though the absolute magnitudes of conduction, 
heating, and cooling rates may differ for H, D, and He, their functional de- 
pendencies and regions of importance are nearly the same in contrast with the 
very different behavior of CO-,. The heat sources used in the computations in- 
clude solar EUV and solar wind heating as well as the important heat sources 
from the electron and ion gases. 

Since C0-, has a very small scale height both the heat input and loss is 
* mainly confined below 150 km altitude with the result that the CO^ exospheric 
temperature is determined almost entirely by the properties of this region. 

Tiie light atomic species absorb energy in a comparitively uniform manner with 
altitude because of their large scale heights above the turbulent region in 
the atmosphere. For these gases, the heat energy must be conducted downward 
into the region of strong coliisional coupling where it can be either trans- 
ferred to the C0 o and lost by infrared radiation or conduction into the 
dense atmosphere below. The details of the heat energy flow depend on whether 
the temperature of the light atomic species is greater or less than the CO., 
temperature. Above approximately 150 km the coliisional coupling is weak so 
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that the neutral gas temperature for He, H, and D is in part determined by 
tiie relative abundance of these species with respect to CO^. An example of 
the temperature separation between the light and CO ^ gases is shown in Figure 2 
for the diurnal variation of temperature and density in the Martian ionosphere 
at 200 km altitude corresponding to the conditions seen by Mariner VI at the 
latitude 4 0 north. The data points for Mariner VI are shown by the small 
rectangles labelled M6 and are clearly in good agreement with the computed 
ionospheric model (see Appendix) . The diurnal neutral gas temperature varia- 
tion is shown by the curves labelled T and T t . For these conditions, T 
(the CO., temperature) becomes larger than T^ (the common H, D, and He tempera- 
ture) shortly after sunrise and stays larger until well after sunset. From 
Figure 2 the time constant for thermal equilibrium is estimated to be about 
three hours. A similar plot for the diurnal variation on Venus differs only 
in the. magnitude of the temperature separation and the lack of noticeable 
asymmetry about local noon because of the very slow rotation rate of Venus. 

The temperature separation shown in Figure 2 has an obvious bearing on 
the problems of planetary evolution through the escape rate of hydrogen and 
on the present altitude distribution of hydrogen. For example, the effect 
of escaping hydrogen is illustrated in the relative behavior of H and D 
shortly after sunrise. On Mars with its weak gravity the normal exobase tem- 
perature is high enough so that hydrogen is depleted and D + is able to dominate 
H + throughout the day. The addition of 0 + to the ionosphere does not change 
their relative behavior. The rapid loss of H + and D + and the low density of 

3 3 

He + (less than 10 an” ) as the day progresses is due to the mutual charge 
exchange reactions between the various ions and their neutral counterparts. 

The net result being that at 200 km altitude most of the produced ionization 
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gocs into the creation of additional CO^ . At higher altitudes other ions 

can dominate (1I + , D + , or 0 + ) except for He + which remains a minor ion by 

charge exchange with all other neutral gases in the atmosphere. 

On Venus the temperature separation can be critical for hydrogen escape. 

Even a 50 °K change in the gas temperature would be sufficient to substantially 

deplete the exosphere of Venus of hydrogen gas. An evaluation of the Jeans - 

1/ 2 

effusion velocity (bulk velocity) V. = (kT L /2T^jg) (1 + R/H) exp (-R/H) 
shows a strong temperature dependence (Vj - 3.017 (T/500) ) in the tem- 

perature range between 500 and 700 °K. R = the radial distance to the exobase, 
ii = kT^/mj.g, m^ = the mass of hydrogen, and g = the acceleration of gravity. 
Possible daytime temperature profiles for Venus and Mars are presented in 
Figure 3 showing a substantial temperature separation in the Venus exosphere. 
The magnitudes of the temperatures and their separation are controlled by the 
.mount of light njutral gases present in upper atmosphere (above about 150 km). 
oir.ee He + and He are free parameters at the present state of knowledge, 
the concentration of neutral helium below 100 km was varied to investigate 
its influence on the temperature and escape of atomic hydrogen, Figure 4. 

As the proportion of helium in the atmosphere below 100 km is increased 
the amount of solar EUV energy available for heating increases for the light 
neutrals and decreases for CO^- In response to this, the light neutral 
temperature T n (H) rises while that of CO 2 decreases. The increase in T n (H) 
at 150 km is limited by the collisional coupling to CO^. At first T^OO 
increases rapidly but as T n (H) approaches ^(CC^) and then exceeds it, the 
heat transfer from CO ^ decreases and finally becomes a heat loss. At 500 km, 


* A least squares fit to V^ with a correlation coefficient greater than 0.99 
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however, where the collisional coupling is very small, the T n (H) and T n (C0 2 ) 
are uncoupled and T n (H) increases exponentially as a function of the atmospheric 
helium content. The effect on the hydrogen density N(H) is shown in Figure 5 
as a power law decrease due to the increased hydrogen escape flux.* The 
derails of the behavior of T (H) and N(H) as shown in Figures 4 and 5 depends 
strongly on the degree of turbulent diffusion in the atmosphere below the 
turbopause region. As the degree of turbulent diffusion increases, as measured 
by an eddy diffusion coefficient (K^^) , the minor constituents are forced 
to follow the C0-, scale height to greater altitudes, McElroy and Hunter. (1363) 
and Herman et. al (1971). This results in a decrease in the density of light 
neutrals above the turbopause and a consequent decrease in T n (H) and an in- 
crease in T^CCO^), Figure 6. Furthermore, when is small T^COp > T^Qi) 

at 150 km and T n (CC> 2 ) < T (H) at 500 km with a relatively small temperature 
separation, but for large T n (H) decreases sharply and the temperature 

separation becomes very significant. The corresponding changes in the hydrogen 
densities are shown in Figure 7. As expected the decreasing T^(H) as 
increases leads to a decrease in the hydrogen escape flux and to an increase 
in \ T (H). Further increases in K e jjy, cause a decrease in N(H) because the de- 
creasing outward bulk flow of hydrogen is eventually dominated by the turbulent 
mixing (i.e. the increase in altitude of the turbopause region). 

For the case of the earth a similar temperature separation occurs between 
the H and He and the other atmospheric constituents (e.g. N 2 , 0 2 , and 0). How- 
ever, the separation turns out to be small since the principal constituent for 
heat loss by infrared radiation, atomic oxygen, extends to great heights and 
provides good collisional coupling with H and He well into the region v.nore 
the neutral gas temperatures are nearly isothermal. 

- 0.22 


* [H] soo = 2.1 x 10 [HeJ 


100 
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In conclusion, it has been shown that a multiple temperature neutral 
atmosphere must be considered for Venus and Mars (and to a much lesser extent 
on Earth) in the interpretation of data obtained from planetary probes. The 
temperature differences between the light atomic species (II, D, and lie) and 
CO, are significant and, particularly on Venus, are very important in deter- 
mining the escape flux of atomic hydrogen and its density above the turbopause 
region. In addition, the Lyman a determination of hydrogen density and tem- 
perature at great heights, Barth (1968) leads to an estimate of the Venus 

s _ 3 

helium concentration at 100 km of 5 x 10 cm and a small value of the eddy 

5 2 .1 

diffusion coefficient (less than S x 10 cm sec ) . 
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APPENDIX 


The computations performed in this work were obtained by simultaneously 
solving the equations of heat transport for each of the electron, ion, and 
neutral gases present in the atmosphere in the manner described by Herman ei. al 
(1971) and Herman and Chandra (1969). In addition the time dependent momentum 
and continuity equations were included to obtain the diurnal variation of 
densities and temperatures presented in Figure 2. A brief summary of these 
equations is as follows. 

For each species, i, the momentum and continuity equations have the 

form 


3V. 3V. 

M. N. + V- 
l l 3t l 3z 


8P. 

- N.M.g + e.N.E - IK.. (V. - V.) 
3z l i b li . ij ' i y 


3N. 

] 

3t 


dZ 


iV.N.) = R. - L d . 
v l i' l Pi 


where the symbols all have their standard meaning. The quantity Kvj can be 
written in terms of collision frequency v— as K-j = v^j where u — 

is the reduced mass. The electric field E is determined from the momentum 
equation for electrons (i.e. -*■ and e^ -*•' |e |). The ion production 

function is the sum of ionization sources from the solar EUV, solar wind 
proton penetration of the ionosphere directly and via charge exchange with 
the ambient neutral atmosphere at great heights, and the mutual charge exchange 

between the various ions and the atmospheric gases. The loss function 

is the sum of the appropriate charge transfer and electron recombination terms 


for each ionic species. 
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FIGURE CAPTIONS 


1 - Venus Dayside and Nightsidc Ionization Profiles. Taken from Fjeldbo 

and Eshlcman (1969) 


2 - The diurnal variation of temperatures and densities for Mars cc puted 

for conditions corresponding to Mariner VI. (4° N latitude, solar 
zeni angle 57°, and orbital longitude 2S4° with solar fluxes scaled 
to • >rrespond to 7/31/69) . The time scale is subdivided into 20 Mars 
hours of 4431 seconds each so as to divide evenly into a Mars day of 
SS620 seconds. On this scale hour 13 is equivalent to an Earth hour 
of 15.6. Solid rectangles labelled M6 represent Mariner VI data points 

5 - Daytime neutral gas temperatures for H, D, He, and CO^ for Venus and 

Mars. Solar zenith angles and solar activity are the same for both 
planets. 

4 - Daytime Venus CC 0 and H temperatures as a function of the helium 

concentration iiVthe lower atmosphere (below 100 km) . The dashed 
lines correspond to an altitude of 150 km and the solid lines to 
500 km. 

5 - Daytime Venus H densities at 150 km (dashed lines) and 500 km (solid 

lines) as a function of neutral helium concentration at 100 km altitude 

6 - Daytime Venus C0 ? and H temperatures as a function of the eddy 

diffusion coefficient. Dashed lines are at 150 km and solid lines 
at 500 km. 

7 - Daytime Venus H densities at 150 km (dashed lines) and 500 km (solid 

lines) as a function of the eddy diffusion coefficient. 
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